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Abstract 

New integrable variant of the one-dimensional Hubbard model with variable-range corre- 
lated hopping is studied. The Hamiltonian is constructed by applying the quantum inverse 
scattering method on the infinite interval at zero density to the one-parameter deformation 
of the L-matrix of the Hubbard model. By construction, this model has Y(su(2))©Y(su(2)) 
symmetry in the infinite chain limit. Multiparticle eigenstates of the model are investigated 
through this method. 
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1 Introduction 



As a model of strongly correlated electrons, the Hubbard model has been attracting much 
interest in solid state physics. Especially, in one dimension, the model is exactly solvable 
and its thermodynamic properties can be calculated out, which give a good testing ground for 
theories of strongly correlated electron systems. From the point of view of the integrability of 
the one-dimensional Hubbard model, there have been many works, including the pioneering 
work of the coordinate Bethe ansatz by Lieb and Wu |p, the quantum inverse scattering 
method g, |, |, g, its S0(4) invariance [0, |, |], Y(su(2))eY(su(2)) invariance in the infinite 
chain limit [l^, and the recent development of the algebraic and analytic Bethe ansatz [11, 12]. 



One of the novel properties of the i?-matrix R{X, fi) associated with the one-dimensional 
Hubbard model is that it is thought to be impossible to express it as a function of a difference 
of two spectral parameters A and /U. This lack of the "difference property" have prevented 
us from investigating underlying integrable structures of the model. For example, it is not 
known whether this i?-matrix is expressible as an intertwiner of a certain algebra. Since the 
methods for calculating various correlation functions known so far [13, 14] requires under- 



standing of such underlying structures of the model to some extent, it is necessary to deepen 
our knowledge of the mathematical structures of the Hubbard model in order to calculate 
correlation functions. 

The absence of the "difference property" is not a mere harm for us. As was noticed by 
Shiroishi and Wadati, it produces a one-parameter integrable extention of the L-matrix and 
the Hamiltonian of the Hubbard model [^]. This extended Hubbard model is interpreted as 
an electronic model with on-site and neighboring-site interactions and correlated hopping to 
the neighboring sites. Though the form of the Hamiltonian is complicated and thus is difficult 
to be equipped with physical meaning, it can be of some help us with the understanding of 
the structures of the original Hubbard model. 

On the other hand, we have recently discovered that the R- and L-matrices of the one- 
dimensional Hubbard model can be put into a formulation of quantum inverse scattering 
method (QISM) on an infinite interval |l^, [l^, which has been applied to other integrable 
model [^, |l^. Through that method, we can derive the existence of Yangian symmetry 
Y(su(2))©Y(su(2)) and construct n-particle states upon zero-density vacuum. Based on 
this work, the aim of the present paper is to put the one-parameter deformed L-matrix, 
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which is described in the previous paragraph, into the same formulation of the QISM on an 
infinite interval. Through this procedure, a new electronic Hamiltonian with variable-range 
correlated hopping arises. It can be embedded in a family of an infinite number of commuting 
operators and thus is interpreted as a one-parameter integrable deformation of the Hubbard 
chain. As is the case for the usual Hubbard chain, Yangian invariance of the Hamiltonian 
and construction of multiparticle states can be directly established as a byproduct of this 
method. 

This paper is organized as follows. In Section 2 we shall explain integrability of the 
Hubbard chain of finite length and a one-parameter deformation of that model. Section 3 is 
devoted to the application of the QISM on an infinite interval to the one-parameter defor- 
mation of the L-matrix. Its resulting new Hamiltonian and commuting conserved operators 
are developed in Section 4. The Yangian Y(su(2))©Y(su(2)) invariance of the model follows 
by construction. In Section 5, we shall construct multiparticle states upon the zero-density 
vacuum by use of symmetries and algebras of some operators. Section 6 contains concluding 
remarks and discussions. 



2 Hamiltonian and Monodromy Matrix on the Finite Interval 

The Hamiltonian for the one-dimensional Hubbard model is 

^ = - E + + ^E (^iT - \) [^ji - 1) > (2-1) 

where cj^ and cj^ are respectively the fermion annihilation and creation operators which 
satisfy the usual anticommutation relations, and njo- = Cj^cja is the particle number operator. 
This Hamiltonian has a large symmetry. First, it is invariant under the partial particle-hole 
transformation: 

c.-T^c,t, c,^^{-iyc]^, U^-U. (2.2) 

Second, it is invariant under S0(4) algebra generated by 5" and r/° (a = x,y,z). S"^ is the 
usual spin-SU(2) operator defined by 

S'' = \E</jaCjp^ (2-3) 

3 

and r]"" is obtained by performing the partial particle-hole transformation (|2.2|) to 5". If we 



go to the infinite chain limit, the symmetry is enhanced to Y(su(2))©Y(su(2)) Yangian |1C]. 



4 



The integrability of the fermionic Hubbard model ( |2.1| ) is based on a local exchange 
relation |^] 

n{x,fi)[Cj{x)(g)sCj{fi)] = [Cj{^x)(^sJOj{X)]n{x,^l). (2.4) 

where (8)s denotes the Grassman tensor product 

[A 0s BU^^ps = (2.5) 

with the grading P{1) = -P(4) = 0, P{2) = P{3) = 1. The expressions for the matrices TZ 
and £ are presented in Appendix A .*| In these expressions, we use a function h{X) defined by 

smh2/i(A) U 

sin2A = 4 • ^'-^^ 

The Hamiltonian (|2.1|) is reproduced by a logarithmic derivative of the transfer matrix TmniX); 



H = -^lnTmn{X)\x=o, T^n(A) = str(r„„(A)) =tr{{a' 0a')%nn{X)), (2.7) 
where the monodromy matrix T is given by 

%nn{X) = Cm-l{X)Cm~2{X) ■ ■ ■ Cn{X) {m > n). (2.8) 

There is known to be an integrable spin chain equivalent to the Hubbard model [P, ^ . If 
we apply the Jordan- Wigner transformation 

Cj'f = a^---aj_^a~, c^i = « ■ ■ ■ CT^_2cr^_iX • • • r/.^r^", (2.9) 

where a and r are the Pauli matrices, and cr^ = (crj it ia^)/2, = {tj ± iTj)/2, we get an 
equivalent spin model 

H = TMt^^J + + YMt^^i + vi) + 7 E -I-/- (2-10) 

j 3 j 

Its integrability is supported by the spin-chain counterpart of the exchange relation 

R{X,^I)[Lj{X)(^Lj{^x)] = [Lj{fi)(g)Lj{X)]R{X,fi). (2.11) 

By using i?i2(A, fi) = PuRuiX, /i), where P12 is the transposition P{x (8) y) = y x, it is also 
written as 

Ri2iX,fi)LijiX)L2j{fi) = L2jifi)Lij{X)Ri2iX,f^), (2.12) 

*The matrix TZ in Ref. |l^, is written as TZ in this paper, following the standard notation. It should 
not be confused with the matrix written in the standard notation as 7?. = VIZ. 
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As Olmedilla et al. Q found, the exchange relation of the fermionic model, ( |2.4| ), and that 
of the spin chain, pH] ), can be transformed into each other. 



One of the peculiarity on the integrability of the Hubbard model known for years is that 
the Tt-matrix, or equivalently the ^-matrix, is believed to lack the difference property, i.e. it 
is not a function of A — /i, nor can it be expressed as /(A) — /(/x) with some function /. This 
has been an obstruction for further investigations of the underlying mathematical structures 
of the Hubbard model. But on the other hand, the lack of the difference property allows 
us to consider a one-parameter integrable deformation of the Hubbard model, as noticed by 



Shiroishi and Wadati |15]. This works as follows. Since the i?-matrix is shown to satisfy the 



Yang-Baxter equation of the form 

i?12(A, fi)Rl3iX, Z^)i?23(^, l^) = i?23(/i, 1^)^13(A, J^)i?12(A, (2.13) 

a new L-matrix defined by 

ii(A). = i2i3(A,z.) (2.14) 

satisfies an exchange relation 

i?i2(A,^)Li(A),L2(/x). = L2(M).ii(A),i2i2(A,/i). (2.15) 
By using Ri2{X, fj.) = Pi2^i2(A,^), it can be alternatively written as 

R{X,fi)[L{X)^0L{fi)^] = [L{fi)^0L{X)^]R{X,fi). (2.16) 

Considering that L{X)u=o oc L{X), which can be checked by a direct calculation, we can say 

that L(A) 1/ is a one-parameter deformation of the L-matrix L[X) of the original Hubbard 



model. This new L-matrix can be used to produce a new Hamiltonian [15|. By using the 
monodromy matrix given by 

TmniX)u = Lm-liX)uLm-2iX)u--- LniX)u {m > u), (2.17) 

the new Hamiltonian is given by 

Hu= ^lntr(Tmn{X)u) 

dX x=u 

^ 4cosh2/i(zy) ^ (^'"■'■^ """^^ " ~ ^"^^^^ ^'""^ ^ ^'"^ 2z.(c]^Cj+i| - c]^i^Cj|) - cos 



(^riji cos^ u — 2nj+ij sin^ u + sin 2i/(c^j^Cj_|_ij — c^j_^-^^Cji) — cos . (2.18) 
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We have performed the Jordan- Wigner transformation to get a fermionic Hamiltonian. The 
reason for choosing the special value A = is to obtain a local Hamiltonian. This Hamiltonian 
is Hermitian if v is pure imaginary, and in that case it represents a model with on-site and 
neighoboring-site interaction and correlated hopping to the neighboring sites. The exchange 
relation ( 2.1(^ ) results in 

[lntrT^„(A)^,lntrT^„(/i)^] = (2.19) 

Therefore, the Hamiltonian is embedded in a family of infinite number of commuting oper- 
ators, and so the model ( |2.18D is integrable. The model ( p.l8| ) includes the Hubbard model 
(12. 1|) as the = case. 



So far we have discussed on a spin-chain version of the exchange relation ( p.l6| ). We can 
employ the same procedure used in Ref. Q to fermionize the exchange relation of the spin 
chain ( |2.16| ). The Jordan- Wigner transformation relates the spin operators and the fermion 
operators as 



m ^' r ' \=n\ ^m- (2-20) 

Here the matrices Vj^ and Vj\^ are defined by 




V,a=\ ' ' ■ (2.21) 




i-i 

k=n k=n k=n 

Then we obtain the fermionic exchange relation; 

lZ{X,lj)[£j{X),^,Cj{fi),] = [£j{fi),^,Cj{X),]lZ{X,fi), (2.23) 

where the fermionic matrix Cj{X) can be calculated from the spin-chain matrix Lj(X) as 

£,(A), = (F,+i,t ® V,+i,OL,iXUVr^' ^ Vf^'), (2.24) 

and the corresponding 7^-matrix is 

7^(A, n) = WR{X, n)W~^, (2.25) 

with W = diag(l, — i, — i, 1) (8> l2- Here denotes the n x n unit matrix. The explicit 
form of Cj{X)i, is presented in [Appendix A . The exchange relation ( p. 16 ) and ( p. 23 ) can 
be considered as a one-parameter deformation of (2.11) and ( |2.4D , respectively. Note that 
R{X, n) or TZ{X, /x) is unchanged by this one-parameter deformation. 
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3 Passage to the Infinite Interval 

Next we pass to the infinite interval limit using the new exchange relation ( p.23 ). The method 



is identical with the one in our previous works on the original Hubbard model 17]. Let 

%nnWu = ^m-l{^)u^m~2i^)iy ' " ^nWu, (3-1) 
%^^{X,^l)i, = TmnWu ®s'Tmn{fJ')u- (3.2) 

%nnWu is a monodromy matrix on the finite interval. To consider the infinite-chain limit of 
the monodromy matrix, we should split off the asymptotics of its vacuum expectation value 
for m, —n oo. Hence, this procedure is restricted to uncorrelated vacua, with which one 
can calculate vacuum expectation value of the monodromy matrix. Among four uncorrelated 
vacua, in which the electron density of each spin is either zero or unity, we take the zero- 
density vacuum [0) to renormalize the monodromy matrix. We use the following two matrices 

v{x), = {o\q{xuo), v^^\x,f,), = {o\q{x),^,Cj{fiUo). (3.3) 

(2) 

in order to normalize %nnW and TmniX^fi)^ respectively; 

t(A), = hm V{X)-'^%nn{X)uV{X):, (3.4) 
t(2)(A,/x), = hm v'^'\X,f,)-"'T(^;i{X,f,)M'\X,f,ru. (3.5) 

m,— n^oo 

(2) 

These limits converge in the weak sense, though the matrices TmnWv and Tmn{X, do not 
have a definite limit when m, —n oo. We shall call T[X)i, a monodromy matrix on the 
infinite interval. It allows an alternative definition: 

f{X), = h + - h) + E(4- - h){CiWu -h) + ---, (3.6) 

where 4(A), = V{X)-^-^q{X),V{X)l. 

For practical calculations, one should be careful that ^'•^-'(A, fi)iy is not equal to the tensor 
product V{X)u (Xis ^{tJ')u- There appear additional off-diagonal elements due to normal order- 
ing of operators. Direct calculations lead us to the resulting forms for V{X)u and V^'^\X, 

V{X)^ = diag{-p8{X,u), pg{X,u), pg{X,iy), -pi{X,u)). 

As for the matrix V^'^\X, p)^, its diagonal consists of the elements of ^(A);^ 0s and 
its non-vanishing off-diagonal elements are 
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V^^\X,p)l''^^' = -psiX,,.)psip,u). 

Note that ^''^•'(A, p)i, is upper triangular. Since the diagonals of V^'^\X, p)y and V{X)®sV{p)v 
are identical, y(2)(A,/^) 

jj can be diagonalized by an upper triangular matrix C/(A,/i) whose 
diagonal elements are all unity; 



y(2)(A,^), = u{x,p),{v{x), ®s v{p),)u{x,p)-\ 



(3.7) 



Direct calculation leads us to a remarkable and surprising fact; t/(A, /u),^ is independent of v. 
It is equal to U (A, p)u=o = U (A, ^u), which has appeared in the analysis of the usual Hubbard 
chain [16|, so we will hereafter suppress the subscript u in U{X,p)u. Its matrix elements are 



U{X,p)i2^2i = U{X,p)i3^3i = -i/92/pio, U{X,p)u^23 = "^(A, /i)i4,32 = ipe/ps, 

U{X,p)24.,42 = ^^(A,/U)34,43 = ip2/ P9, U{X, p) 23,41 = -U{X,p)z2,Al = ipe/P?, 
C/(A,^)l4,41 = -p^/pi, 

where pi = Pi{X,p). 

Taking the vacuum expectation value of the local exchange relation (2.4) yields 

n{x,p)v^^\x,p), = v''^\p,x),nix,p), (3. 



and we conclude that 



7^(A, /i)t(2) (A, p), = t(2) (p, A),7^(A, p) . 



(3.9) 



Finally, collecting ( |3.9D and other equations together, we arrive at the exchange relation 
for the monodromy matrix T{X)y on the infinite interval, 



where 



n^+\x,p) \f{x\®,f{p)J\ = [t(/i), 0s t(A) J n^-\x,p) 



n^^\x,p), = u±{p,x)-'n{x,p)u±{x,p). 



U±{X,p),= hrn v'^^\X,p)-"'[V{Xr(^sV{p) 

m— >±oo 



(3.10) 

(3.11) 
(3.12) 
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Since the calculation of the matrices U±{X, and '^'•^^(A, is rather technical, we do not 
reproduce it here. Its details are presented in Ref. |17] in the case of = (Hubbard model). 
The only point we should note here is that apart from some singular points (e.g. X = fi), we 
can say that U{X,fi) = U±{X,^)u and 

1z{x,^l) = 1z'•^\x,^),= 
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Pi) 



(3.13) 



(3.14) 



where pi = pi{X,fi). Let us write the elements of T(A) as 

( Dn{X)u Cn(A), CuiX), DuiX), \ 
Bu{X), An(A), Ai2(A), BuiX), 

B2l{X), A2l(A), A22{X), B22{X), 
\ D2l{X), C2l{X), C22{X)u D22{X), ) 

Since the 7^-matrix is independent of the value of the commutation rules between the 
elements of 7"(A), which are obtained from the exchange relation ( |3.10| ), are completely the 
same as the ones in the u = case, i.e. the usual Hubbard model. The complete list of the 



commutation rules are found in Appendix B of Ref. |17|, and it is the same in the present 
case. 
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4 Yangian Symmetry and Commuting Operators 

If we follow the notion of the quantum inverse scattering method, the remaining task is to 
investigate the meaning of each matrix element of T{X)y. As is also the case for the Hubbard 



model [jl^l, the commutation relations between the elements of the submatrix A{\)y decouple 



from the rest of the algebra. 

r(A, li) {A{\\ A{ii),) = {A{ii), ® A{\),) r(A, (4.1) 

where 

^ PSPi - P2 + P9P10'P ( (\ \\ (An\ 

r{X,n) = [pj = Pj{\p)), (4.2) 

P4iP3 - Pi) 



and P is a 4 X 4 permutation matrix [Vx ® y = y ® x). As is remarked previously, (4.1) is 
identical with the one in the Hubbard model {u = 0) and we can follow the same argument 
as in the previous work |17|. By the reparameterization 



v{\) = -2i cot 2A cosh 2/i(A), (4.3) 
the iZ-matrix r(A, p) turns into the rational i?-matrix of the XXX spin chain, 

nA.M) = '^,r'°'''"'"f'r . (4.4) 



Let us expand A{X)y in terms of v{X) 



-1 



oo -, / 3 \ 

mv = /2 + if/ E -^7X^1 E Qn(^)^'^ + ^?n(^)^2 
n=0 ^ ' \a=l / 



(4.5) 



where = — o"^, fj^ = o"^, and = a^. This expansion can be achieved by considering the 
limit v{\) — > cxo as 9(A) — > oo and by choosing the proper branch of solution of eq. (2.6), 
which determines /i as a function of A. (|2.6|) implies that 



g-2MA) = _^ sin 2A ± yi + J sin . (4.6) 
To achieve convergence of the matrix elements A^X),^, we have to choose the lower sign 



here. Then it follows from general considerations pQ, 22] that the first six operators 



Qo(^)i Qii^) generate a representation of the Y(su(2)) Yangian quantum group. 



There is an alternative description of the Yangian Y(su(2)) | p4| described below. The 
Yangian Y(su(2)) is a Hopf algebra spanned by six generators Qn (^ = 0)1) = x^y^z), 
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satisfying the following relations, 

VO' HQ — J WQi 



Qi,Qi , Qo,Qi 



+ 



Qi,Qi , Qo,Qi 



^2^^abkefgjcdk ^ ^cdkefg jabk^i^QC^ g/^ 



(4.7) 
(4.8) 

(4.9) 



Here k is a nonzero constant, /"^'^ = ie"'^'^ is the antisymmetric tensor of structure constants 
of su(2), and A''^'"^''^ = f^dkjbeijcfmjkim^ rpj^g bracket { } in (|^ denotes the symmetrized 
product 



{xi,X2,X3} = X! XaiX^2Xa3- 



(4.10) 



The Hopf algebra structure of Y(su(2)) is described in ref. |24] and its representation theory, 
which will be used later, is developed in Ref. |2^, 26 1. 

We can use ( [4 .51 ) and ( p.6| ) to get the representation of Yangian generators; 

IT.<A»'^J,P' (4-11) 



2 sin cos u 



^^(tan i.)-^e2^('^)(2-"-"^)<^cLc,^ 



i>j 



+ cot ^)^-^e-2'^('^)(^-"-"^)<^cLc,^ 



i<j 



iU 



— Y,sgn{j -i)a^pcl^cl^Ci,^Cj^f3 + iUcot2h{i^) sin^ Qg(i^). (4.12) 



In this case the constant k in (|4.9| ) is equal to iU. Note that Qol'^) = "S*" is just the operator 
of the a-component of the total spin. The representation of the Yangian algebra in the usual 
Hubbard model |1C, 16| is a special case of = in ( 4.11| ) and ( [4.12| ). 



Since the quantum determinant 

DetgA(A)^ = ^ii(A)^A22(A)^ - ^i2(A)^^2i(A)^, 
where f (A) = f (A) — iU, is in the center of the Yangian 

[BetgA{X)^,A{fi)^] = 0, 
it provides a generating function of mutually commuting operators, 

[BetgA{X)^„BetgA{n)^] = 0. 



(4.13) 



(4.14) 



(4.15) 



12 



The asymptotic expansion in terms of v{X) ^, 



Det,A(A). = l + if/|:-^, 



produces Jo{i') = 0, Ji(z^) = i^^iong, where 



long 



sm cos ly 



i>j 



anz^)*-^e2'^(^)(i-"»---"^-.-)4cj<. 



cot ;,)»-ie-2'»H(i-"'.--"j-.-)cLc. 



^7(1 + 2sin2 



riji - - I I n, 



1\ 1' 



(4.16) 



(4.17) 



Due to the relation ( [4.15 ), the J„(zv)'s mutually commute. Therefore, H\ong can be em- 
bedded in a family of infinite number of commuting operators, and can be regarded as a 
integrable Hamiltonian. Moreover, (|1|) indicates the Y(su(2)) invariance of the model; 



[Qg(i/), Hionsl = = [Q?(z.), i^iong], (a = 1, 2, 3). 



(4.18) 



Especially, it implies that the model is su(2) invariant. By subtracting a constant from i^iongj 
we can make this Hamiltonian invariant under partial particle-hole transformation (p.2|): 



Hi 



long 



sm v cos v 



^(tanzy)'^-^e2'^('')(^-"''--"^'-)cLc. 



i>j 



cot z.)i-ie-2'»('^)(i-">---"i.-)cLc 



+ [/(l + 2sin2 i/)^ 



n 



(4.19) 



This complicated Hamiltonian can be made simpler by noting (p.6[) to get 



H" 

long 



-i sin u cos i^Hi 



' 5Z ( ». 



long 



+ 



6 - 2r2 



sm 



2; ' 



(4.20) 



where r = icot u, J = —2\h{u) and r and J are real. We will, however, use the Hamiltonian 
( 4.19| ) instead of ( 4.201 ), since it is easier to extract informations of the model out of the 
quantum inverse scattering method. 
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The Y(su(2)) invariance of H^^^^ shown in ( [4.18 ), together with the invariance under the 



partial particle-hole transformation (2.2), leads us to the result 



[Qg(i.),i?^g = = [QU'^),HU^], (4.21) 
[QS'(z.),i7Ug] = = [Qrii^),HU,l (4.22) 



where Qn{^) is obtained by performing the partial particle-hole transformation (2^) to 
Qn{^)- By construction, the operators Qni^) form another Y(su(2)) algebra, and we can 
straightforwardly verify that [Qm{^)^Q^n{^)\ vanishes for a,b = x,y,z; m,n = 0, 1. Putting 
all things together, we can say that the Hamiltonian H{ong Y(su(2))©Y(su(2)) invariant. 

Hereafter we shall assume -ff/ong to be Hermitian, i.e. U is real and both u and h{h') 
are pure imaginary. With this assumption, H[^^^ can be regarded as a new Hamiltonian for 
electrons with on-site interaction and variable range hopping. The amplitude of the hopping 
decays exponentially with the hopping range. There is also an interference effect due to the 
term exp(ib2/i(i/)(l — nj^_o- — nj^_CT)). One can easily check that in the u = limit the hopping 
terms vanish except for the ones to the neighboring sites, and the usual Hubbard model is 



restored in this limit. In that sense it is an integrable extension of the Hubbard model (2.1 
with variable range hopping. 

There is another integrable Hubbard model with variable range hopping discovered ear- 



lier [27|. Its Hamiltonian is given by 



H=J2 tii-j)4.c,a + Uj2n^^n,^, (4.23) 

j 



with 



t{s) = -it(-l)" - sin — (4.24) 
vr L 



or 

t{s) = -isinhK(-l)7sinh(Ks). (4.25) 

Though both our model and the above model contains the usual nearest-neighbor-hopping 
Hubbard model as a limiting case, we do not know whether they can be related with each 
other. 

In closing this section, we shall add a comment. The Hamiltonian i?iong obtained here is 
different from Hu obtained from a logarithmic derivative of the monodromy matrix %nn{^)u 
on the finite interval. Such things do not occur in the previously studied cases of the fermionic 
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nonlinear Schrodinger model |22] or the Hubbard model [16, 17]. The relation between the 
two Hamiltonians -ffiong and Hi, is left as a future problem. 



5 Construction of Eigenvectors 



5.1 Creation Operators of Quasiparticles 



As is the case with the usual Hubbard chain and other integrable models, the entries of the 
monodromy matrix T(A) can be used to construct multiparticle eigenstates on the vacuum. 
By calculating commutators between these entries of T(A) and the particle number operator 
N, we can see that -Bai(-^)i/ and C2a{^)u add N by one and D2i{\)u adds N by two, while 
Ba2{^)u and Cia{X)i, reduce N by one and Di2{X)u does by two. Aab{X)u, Dii{X)i, and 
D22iX)u keep N unchanged. To gain further insight, let us calculate actions of the operators 



in T(A) onto the vacuum |0). By using (3.6), some simplest ones are calculated as follows; 



BiiiXUO) 
S2i(A)^|0) 

C2l{X)u\0) 
C22(A).|0) 

D2i{XU0) 



1/06 
P9 



^^e-^^'P(^'^)4^|0), 



P2 
Pi 



-iifc(A,v) t 



Pi 



E "^jT^U 



P9P1 



Pi 



|0>, 



(5.1) 
(5.2) 
(5.3) 
(5.4) 

P9P1 

(5.5) 



where pj = pj{X^v) and 



-'^(^'-) = -p,{X,y)/pi{X,v), e-^'^^'-) = -p^{X,v)/p,{X,y). 



(5.6) 



The commutators between T)eiq{A{p)y) and the various operators in T{X)y are calculated 
from the exchange relation (|3.1C| ) . The resulting commutators are the same as the case of the 
original Hubbard model (i^ = 0), which is summarized in Appendix B.2 of Ref . |l7|] . Hence, 
the commutators between the Hamiltonian ifiong and the operators in the matrix T(A)jy are 



15 



also the same as the u = case; 



[Along, Bai{X)u] = -{2cosp{X) + U/2)Bai{X)u, (5.7) 

[Hiong,Ba2iX)u] = (2 COS A:(A) + C//2)S,2 (A)„ (5.8) 

[Hiong,Cia{XU = (2cosp(A) + C//2)Cia(A),, (5.9) 

[^iong,C2a(A).] = -(2cosA;(A) + ?7/2)C2a(A)., (5.10) 

[Hions,D,2{XU = 2(e'PW+e-''=W)I)i2(A),, (5.11) 

[Hiong,D2i{XU = -2(e^PW+e-^'=W)Z)2i(A)., (5.12) 



where 



e 



Other entries Z?ii(A), I?22(A) and ylab(A) commute with -ffiong- The above results justify 
the interpretation of Bai{X),C2aiX) and -D2i(A) as creation operators. Bai{X) and C2a(A) 
create single particle excitations, whereas -D2i(A) creates a bound state of two particles. For 
example, from -ffiong|0) = and ( ^.10 ) we deduce 



i?longC2ax(Al). • • • C2a„(A„),|0) = " ^(2 COS k{Xj) + U/2) C2a, (Al). • • • C2a„ (A„), |0) . (5.14) 

Similarly, the applications of the operators Bai{X) or mixed products of operators Bai{X) 
and C2a(A) on the vacuum produce eigenstates of the Hamiltonian. 

Let us consider the relation between these creation operators and Yangian Y(su(2)). Since 
the Yangian generators Q^i^) = 0, 1; a = x,y,z) are coefficients of the power expansion 
of A{ij,),^, the commutators between Qnii^) and operators B{X)y, C{X)y, and D{X)y can be 
obtained from (|3.10|) . The results are the same as in the case of the Hubbard model ||l^ ; 

[Qg(z.),i3(A),] = -^a«5(A),, (5.15) 

[Qt{v),B{XU = s\np{X)d''B{X), + ^e'''"'a'>B{X)uQl{v), (5.16) 

[Ql{v),C{XU = ^C(A),a^ (5.17) 

[Q1{v),C{XU = -sinA:(A)C(A),<7'^ + |e"''^C7(A),a^QS(z.), (5.18) 
[Ql{u),D[XU = [Q'i[u),D{XU=Q. (5.19) 

These commutators will be used to investigate Yangian representations of the multiparticle 
eigenstates. 



16 



5.2 Scattering States 



Observing the cases of other integrable models studied ear Her |18, 28], we propose the fol- 



lowing two pairs of normalized creation operators of scattering states, 

R^{X)l = i^-'-Pli^C2a{X)uD22{X)-' {a = l,2), (5.20) 

Ra{X)i = i"-i^4v4^3-a,i(A).I)ii(A);i (a = l,2). (5.21) 

In these formulae a = 1 corresponds to spin- up and a = 2 to spin-down. The numerical 
prefactors have been obtained by demanding that Ra{X)l and Ra{X)l generate normalized 
one-particle states. 



Ra{X)l\0) =J2e-'^'^''''^c]jO), Ra{X)t\0) = 5]e-^-p(^'-)4j0). (5.22) 
j j 



Hereafter we assume that A is chosen in such a way that Ra{X)l, and RaiX)}^ create physical 
states. This means for Ra{X)}^ that k{X, v) has to be real and for Ra{X)\, that p(A, v) has to 
be real. 



By the method in Ref. hermitian conjugation can be performed to the operators 
Ra{X)ij and Ra{X)u, and the resulting normalized annihilation operators are 

Ra{X)u = i^-^ P'^/"';"] ZJn(A);^Ci,3-.(A)., (5.23) 

Ra{X)u = i"-'^4T74 ^22(A);iB„2(A)., (5.24) 

with A' = 7r/2 — A*. The commutation rules between the operators i?a(A)J,, i?a(A)J,, Ra{X)iy, 
Ra{X)y are 

R^{X)IR(,{^)1 = -r{X,f,)^s,apR^{f^)lRsiX)l (5.25) 

Ra{X)i,Rfs{fl)l = -r{^,X)^a,5l3R-y{P')lR5{X)u, (5.26) 

Ra{X)lRf,{fi)l = -r{fi,X)^s,af^R^{p)lR5Wl (5.27) 

Ra{X)uRl3{p)t = -r{X,n)ja,Sf3R'yip)lR5iX)u, (5.28) 

RaiX)lRp{fi)t = -Rp{fi)tRa{X)l (5.29) 
RaiX),Rp{f,)l = -i?^(/i)ti?,(A),. (5.30) 



Here we neglected 5-function contributions from some singular points, e.g., A = /i. 
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The operators Ra{X), i?a(A)^ and Ra{X), Ra{X)^ form a representation of the graded 
Zamolodchikov-Faddeev algebra with S'-matrix r(A,/i). These representations may be identi- 



fied as representations of left and right Zamolodchikov-Faddeev algebra, respectively 

31, 32, |3^. The operators Ra{X)^ and Ra{X)^ are graded as odd, which implies that they are 

creation operators of fermionic quasi-particles. 

We shall present two-particle states generated by Ra{X)^ to know physical meanings of 
Rl{X)„; 

Ri{X)iRi{fi)l\0) = ^4^q\e-'^''^(^'-)e-i"=('^'-)|0), (5.31) 
R2{X)iR2i^i)i\0) = ^4^4^e-'^-^(^'-)e-"'=('^'-)|0), (5.32) 



viX)-v{fM)+iU 

+ 5j-je-'-''^'=(^''^)+'=('^'''»F(A, ^, u)] [0), (5.33) 
R.iX)lR,U\0) = EsVt \0ij > Oe---^^^--^e-'^(^^-) ^^;("J-;(^^\ ^ 



?;(A) - z;(/x) + 1^7 

+ ^^•ie-'-''^'=(^''^)+'=('^'^)>F(A, |0) , (5.34) 



where 



F(A, ^, I/) = ^ (A, ^) f eMA)+MM)-2Mz.) cos a cos ;u + e-MA)-M/.)+2M-) gin A sin fi) . (5.35) 

P4P8 ^ ' 

From these wavefunctions, like in the case of the Hubbard model, we conjecture that the 
n-particle state 

i?ai(Ai)t---i?„„(A„)t|0) (5.36) 

is a normalized in-state if A;(Ai, i/) < • • • < A;(A„, and a normalized out-state if A;(Ai, v') > 
■ ■ ■ > k(Xn,'^)- Here "normalized" means that the magnitude of the incident wave is unity. 
Similarly, for the R operators, we conjecture that the n-particle state 

Rar{Xl)i---RaMl\0) (5.37) 

is a normalized out-state if p(Ai,z/) < • • • < p{Xn-,v) and a normalized in-state if p(Ai,z^) > 
• • • > p{Xn, f). We have not found out its proof yet. 
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The Yangian representation of the multiparticle states can be investigated by the fohowing 
commutators from ( |5.15| )-( |5?T9| ); 

[QSH,i?„(A)t] = ii?^(A)ta^,, (5.38) 

[Q?(z.),i?,(A)t] = -sinA:(A)i?^(A)ta^„ + |e-''-i?^(A)t (5.39) 

[Qt{u),R^{\)l] = \Rp{\)lal^, (5.40) 

[Q?(z.),i?,(A)t] = -sinp(A)i?;3(A)ta^„ - ^e'^'^ Rp{\)la'p^Ql{v) . (5.41) 



These formulae induce an action of the Yangian on n-particle states |34, 17 1. Noting that 
Qo(^)|0) = = (z^)|0), we obtain the action of the Yangian on the n = 1 sector as 

QS(z.)ii«(A)t|0) = ^a^„i?^(A)t[0), (5.42) 
QUi^)RUX)i\0) = -sinA;(A)cT^„i?;3(A)t|0), (5.43) 

which is identified as the fundamental representation VFi(— 2 sin A;(A)). 

The 2"-dimensional representation formed by n-particle states ( ^.36D can be studied by 
the similar manner as Ref. ||2^, |l^, and is identified as the tensor product representation 
M^i(— 2 sin A;(Ai)) • • • (8) PFi(— 2 sin fc(A„)) with co-multiplication A defined by 

A(<3S) = Qo^l + l^Qo, (5.44) 
A(Q?) = 1 + 1 ® + C/e^^^Ql; ® Qg. (5.45) 

This representation is irreducible since A;(Aj)'s are real. Hence, we conclude that all the n- 



particle states of (|5.3£| ) can be constructed by applying the Yangian generators Qnii') to the 
highest weight state 

i?t(Ai)t---i?T(An)l|0), (5.46) 

which is clearly proportional to 

ct(A;(Ai,i/))t...ct(A:(A„,iv))t|0), (5.47) 

with Ca{k)^ = J2j ^~^''^'^]a- ^ operators instead of R, we reach the similar conclusion 

with k{X) replaced by p{\) and the definition of the comultiplication changed to 

A'(QS) = Qg (8) 1 + 1 QS, (5.48) 
A'(Q?) = g^f 1 + 1 - ^7e«''^Q^ ® (5.49) 
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5.3 Bound States 

In order to investigate structures of bound states, let us begin with the two-particle bound 
states. Among two-particle states 

C2a(Al).C26(A2).|0), (5.50) 

we should set (a, b) = (2, 1), (1, 2) in order to obtain bound states, which follows from explicit 
calculation of wavefunctions. In the former case, in which the eigenstate is calculated as 

C22(Ai).C2i(A2).|0) oc ^4^4^ [e{j > l){v{Xi) - t;(A2))e-^^''=(^^''^)e-^"=(^^''^) 

+ eU < l)iv{Xi) - v{\2) + if/)e-y'=(^i''^)e-"'=(^2''^) + e{j < Z)(-i[7)e-i"^(^i''')e-^-'''=(^2''') 
+ 5,;e"'^'W^^''')+'=(^^''^»(^;(Ai) -t;(A2) +iC/)F(Ai,A2,i^)] |0). (5.51) 

The condition for it to be a bound state is 

v{Xi) - f (A2) = -iU, (5.52) 
Qk{Xi,i^) = -9/c(A2, 1^) < 0. (5.53) 

Provided that these conditions hold, it follows that 

C22(Al).C2l(A2). = -C2l(Ai),C22(A2)., (5.54) 

which implies that (a, b) = (2, 1) and (2, 1) cases give the same bound state. To summarize, 
among the two-particle states ( |5.5C1| ), there is only one bound state, which is achieved in the 
case {a,b) = (2,1) with the conditions ( 5.52 ), ( 5^5^ ). We have not, however, succeeded to 



investigate this condition further due to the complicated form of the function k{\,v). 

Due to this complicated form of k(X,i/), we have not found out general forms of multi- 
particle bound states or their creation operators. But if we simply mimic the construction 



of bound state operators of the original Hubbard model in Ref. |17|, we can formally make 
"bound state operators" by 



C2^™'''(Ai, • • • , A2m)i/ — C'22(Al)i/C2l(A2)iyC22(A3)i^C2l(A4)jy • • • C22 (A2m-1 )i/C'21 ( A2mfe55) 
-^22 (Al, • • • 5 X2m)u = D22{Xl)uD22{X2)u ' ' ' -D22 (A2m-1 )i/-D22 ( A2m,)i/ • (5.56) 
R^^™'\Xl, • • • , A2m)I, = C2^'"'*(Al, • • • , A2m)v-D22'"''(Al, • • • , X2m)y'^ ■, (5.57) 
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where 



KX2s) + p(A2s-i) = vr (mod 27r), (5.58) 

smfc(A25-i) = sin/i:(Ai) + '^^^^ ^' (5.59) 

The commutation rules of i?(2™)t with i?^^")^', or Q^i^) are the same as those in the case 



of the usual Hubbard model ((6.58)-(6.60) in Ref. |17|). But the serious problem with this 
operator 

^{2m)t is that we do not know whether it certainly produces physical states. As 
the case of two-particle bound states (m = 1) is already difficult to study, there is little hope 
that we can get deep understanding of multiparticle bound state operators. 



6 Concluding Remarks and Discussion 

In this paper we have introduced new integrable variant of the nearest-neighbor Hubbard 
model with variable range hopping. We have constructed it by the quantum inverse scat- 
tering method on the infinite interval at zero density, using the one-parameter deformation 
of the >C-matrix of the Hubbard model. By construction, together with the knowledge of 
the case of the Hubbard model studied earlier, this Hamiltonian is among an infinite num- 
ber of commuting operators and thus integrable. Moreover, it commutes with operators 
Qni^) ("- = 0,1; a = x,y,z), which form a representation of the Y(su(2)) Yangian. If 
we take the Hamiltonian H^^^^ or H{'^^g instead of -ffiongi it is invariant under the partial 
particle-hole transformation and is Y(su(2))©Y(su(2)) invariant. Normalized creation and 
annihilation operators of quasiparticles are explicitly constructed and are shown to form the 
Zamolodchikov-Faddeev algebra. Multiparticle scattering states are constructed with these 
operators, while bound states still require further study. 

The merits of deriving the Hamiltonian along such line are as follows. First, the existence 
of the Yangian symmetry can be established without any ad hoc methods. Second, the forms 
of multiparticle states upon the zero-density vacuum can be derived without any ansatz. 
They are calculated directly from actions of the elements of the monodromy matrix T(A),y 
on the vacuum. Derivation of multiparticle wavefunctions by making some kind of "Bethe 
ansatz" is rather difficult due to the complicated structure of the Hamiltonian. 

Although it is not easy to interpret the physical meaning of the term of the phase factor 
exp(zb2/i(i^)(l — Ui-a — rij-o)) in the Hamiltonian, it would be interesting to investigate 
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thermodynamic properties of this new integrable model and it will be studied in separate 



papers |35 
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Appendix A Expressions of the 7?.-matrix and the >C-matrix 

In this appendix we present the expressions of the 7t-matrix and the £-matrix. The C- 
matrix [0 is given as 



LAX) = e 



/i(A)(T^(gio-^/2 



with the grading P{1) = 0, P{2) = 1 and functions 

fjaW = (icot A)"^"" sin A, gjaW = (-itanA)"^'" cos A. 



(A.2) 



The 7^-matrix is given as 



7t(A,M) 
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(A.3) 



where pj = pj{X,p) is defined by 

(e' cos A cos // + e~' sin A sin //) p2 , 
(e' sin A sin p + cos A cos p) p2 , 
(— e' cos A sin p + e~' sin A cos p)p2, 

(e' sin A cos p — e~' cos A sin p)p2, 

e' cos A cos p — e~' sin A sin p 

2\ ~ ^2, 

cos^ A — sm p 

— e' sin A sin p + e~' cos A cos p 

cos^ A — sm p 

g-2/i(/i) j,Qg A sin A — e"^'*^'''^ cos p sin p 

cos^ A — sin^ p 

P^-Pi, P8 = P3-Pl, 



Pi = 

P4 = 

P9 = 

PlO = 

P3 = 

P5 = 

P6 = 

P7 = 



with / = h{X) — h[p). The function h{X) is defined by ( |2.6| ). These matrices are identical with 
those in Ref. Q, except for the point that the spectral parameters A and p are shifted by 
7r/4. These matrices satisfy the exchange relation (2.4). Its one-parameter deformed version 



23 



is given by 



2.16|) with the T^-matrix presented above and the new >C-matrix given by 







= Pinj\nji - i/)io(nj| + Uji - 
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where pj = pj{X,v). One can easily check that Cj{\)u=Q oc Cj{X). 
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